Ultracold Fermion Cooling Cycle using Heteronuclear Feshbach Resonances 
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We consider an ideal gas of Bose and Fermi atoms in a harmonic trap, with a Feshbach resonance in 
the interspecies atomic scattering that can lead to formation of fermionic molecules. We map out the 
phase diagram for this three-component mixture in chemical and thermal equilibrium. Considering 
adiabatic association and dissociation of the molecules, we identify a possible cooling cycle, which 
in ideal circumstances can yield an exponential increase of the phase-space density. 



Intro duct ion-Feshbach resonances 1] have become 
a powerful tool for controlling atomic interactions in ul- 
tracold gases 5, 6, 71,13, 0- By tuning an ex- 
ternal magnetic field, the energy of a molecular bound 
state can be adjusted to be on resonance with a scat- 
tering state of two free atoms, thus facilitating external 
control of the atomic interactions. Tuning the molecu- 
lar state below threshold can lead to the production of 
stable ultracold diatomic molecules. This has been ac- 
complished both in bosonic and fermionic gases 0. 
In all of these cases the diatomic molecules are homonu- 
clear bosons. Two recent experiments, at MIT 8j and 
at JILA 9] have observed Feshbach resonances in Bose- 
Fermi atomic mixtures. The molecular species in this 
case is a heteronuclear fermion. 

In this work we construct the phase diagram for an 
ideal trapped gas of bosonic and fermionic atoms, and the 
heteronuclear dimer formed from then i llOl . In contrast 
to previous work on Fermi gases [ill . Il2j . we find that 
adiabatic passage from free atoms to molecules can lead 
to cooling of the system under conditions comparable to 
those of current experiments. Based on this observation 
we propose a thermodynamic cooling cycle, consisting of 
adiabatic interconversion between atoms and molecules 
in conjunction with selective removal of atoms. This cy- 
cle is capable of reducing the absolute temperature of the 
system while increasing the atomic phase-space density. 

Equilibrium theory- We consider an ideal gas mix- 
ture composed of fermionic (/) and bosonic (6) atoms 
and heteronuclear fermionic molecules {hf) trapped in 
an anisotropic harmonic optical trap. The trapping fre- 
quencies for the different components are taken to be: 
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where i = {x^y^z}, a\ is proportional to the atomic 
polarizability, and rrik are the particle masses, with 
k = {f^b^bf}. The expression for the molecular fre- 
quency colj^ in terms of the atomic frequencies is ap- 
propriate when the molecular internuclear separation is 
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large compared to the mean atomic radii. It is possible 
to change the ratio cj^/cj^ over a wide range by vary- 
ing the wavelength of the trapping laser, and Ref. 13] 
gives specific prescriptions for doing this. We assume 
chemical and thermal equilibrium: /i/ + /i5 = /i^/ and 
Tf = Ti) = Ti)f = where fik are the chemical poten- 
tials of the three components, and T is the temperature. 
The energy of the molecular state as measured from the 
atomic dissociation continuum is e^es^ which is a func- 
tion of an applied magnetic field. The results of this 
paper are based on considering the phase diagram of this 
three-component system as a function of T, eres, and the 
atomic populations. 

We employ a semi-classical approximation, such that 
populations of each component are given by 
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where the fugacities Zk are 
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Here is Boltzmann's constant, and J^n{z) and Qn{z) 
are the Fermi-Dirac and Bose-Einstein integrals, respec- 
tively 0, ^3 . This approximation is valid provided 
k^T ^ fuJok^ where Ok = {uo^uo^ix)^)^ is the mean har- 
monic frequency. Nh{T^ii})) represents the number of 
non-condensed Bose atoms; we denote the number of 
Bose condensed atoms by A/'c, and Ni^ = Nc -\- Ni^. 

For a given T, eres, total number of atoms A/", and 
boson-fermion population difference AN the atomic 
chemical potentials are determined by the conditions 

N = iV^(T,/ij)+7V,(T,/i5) + 27V,^(T,/i5/,eres),(6) 
AN = Nf{T,f,f)-N,{T,^,). (7) 

For T above a critical temperature Tc all bosons are 
thermal and A/'c = 0, whereas for T < Tc a Bose-Einstein 
condensate (BEC) emerges: Nc > 0^ jUb = {jUbf = M/), 
and 



iVb(T, Mb = 0) = C(3) 
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FIG. 1: Phase diagram for trapping parameters of the JILA 
system p^. with AN — 0; this representation is universal 
in N subject to validity of the semi-classical approximation. 
Left frame: molecule fraction, r/tf, as a function of Cres and 
T; right frame: condensate fraction, r]c. White lines indicate 
adiabatic trajectories, i.e. curves of constant entropy. 



where ({n) is the Riemann zeta function. 

Phase diagrams- We now discuss the phase diagrams 
resulting from the numerical solutions to Eqs. © and 
Q. Defining the molecular and condensate fraction as 
775/ = 2Ni)f/N and r]c = J^c/J^b^ respectively, we plot 
these quantities in Fig. ^ versus temperature and e^es 
for a system where AN = 0, using trapping parameters 
for the JILA experiment 15]. The T and eres axes are 
normalized by T-p/rij and /ceTp/r^/, respectively, where 
Tp is the Fermi temperature of the fermionic atoms, 
A:bTf = hOf{6xNy/^, x = (TV/ + Nbf)/N is the total 
fermionic fraction, r^f = Of/Obf and = Of/Ob. For 
a finite population difference the maximum number of 
molecules is TV^^^ = {N -\AN\)/2. 

When AA^ < 0, there will always be free bosons in the 
system, and condensation is possible for any value of eres- 

(c) 

When AA^ > 0, there is a critical resonance energy, eres, 
at which the condensate fraction vanishes at T = 0. This 
can be seen in Fig. ^ Setting A^5 = in Eqs. © and {TJ 

and = eres + {6Nbf)'^^^hObf at T = we find that 
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In a similar manner we can solve for the critical tem- 
perature Tc for BEC in the limit eres/^B^F ^ 1, where 
Nbf =0. From Eq. © we find 



^ = [6C(3)]-^/^- 
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In the case of equal populations x = 1/2, the limits take 
a simple form e^s = —kBTF/Hf and Tc ^ O.SlSTp/r^. 
These two points define the left and top boundaries of 
the condensate region in Fig. Q 

Adiabatic sweep- During an adiabatic sweep, the 
resonant energy is varied on a time scale much longer 
than the average relaxation time of the system ^] . Then 



entropy is conserved, and the system follows an adiabatic 
path through the phase diagram. The total entropy of 
the system, S{T)^ is given by the sum of the entropies of 
each of the three components: 



Sf{T,^^f) = kBNf 
Sb{T,iib) = ksNb 

Sbf{T, flbf,€res) = k^Nbf 
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The final temperature and chemical potentials at the 
end of the sweep are determined from the condition 
S{T^) = S{T^) together with the constraints Q and JTj). 
Figure n shows the contours of constant entropy in the 
phase diagrams for the case of equal populations. Start- 
ing with a pure atomic system at eres ^ ^b^f we observe 
that there is a temperature, Tiso, such that for > Tiso 
the final temperature after complete association of the 
gas is higher than the initial temperature. However, for 
< Tiso we find the interesting result that the temper- 
ature of the system decreases after the atoms have been 
converted into molecules. This can be understood from 
the fact that at low temperatures, the majority of the 
free bosons are condensed and do not contribute to S. 
By sweeping the value of eres across the resonance we go 
from an almost perfect Fermi sea of atoms to an almost 
perfect Fermi sea of molecules. But since Obf < 00 the 
final temperature must be lower if S{T^) = S{T^). 

An analytic expression for the relation between the 
initial and final temperatures can be obtained in high- 
and low-T limits. For T ^ (Tf,Tc), the entropy of each 
component is of the Boltzmann form [l^ 



Sk{T) = ksNk U - In 
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One obtains T-^ from S{TJ^) = S{T'); in the case x = 1/2 
the result is 



T'\ 



The bosonic entropy for T <C (Tf,Tc) is: 
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while for the fermionic atoms the low temperature ex- 
pansion of Tn{z) |1^ yields: 



S{T) = TT^kBNf 
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The entropy of the molecules is given by a similar expres- 
sion. The adiabatic constraint S{T^) = S{T^) in general 
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FIG. 2: Final temperature after adiabatic sweep from atomic 
to molecular configuration as a function of the initial temper- 
ature (dots). The frequency ratios correspond to those of the 
JILA experiment 15], and x — 1/2. The analytic relation 
(|2()|) is shown for comparison (solid line), along with the line 
rpf _ rp% (^jQ^g]^g(j line). The inset shows the same numerical 
results compared with the analytic expression (|15|) in the high 
temperature limit (solid line). 



FIG. 3: Ratio /T' after adiabatic conversion of atoms to 
molecules (Eq. (fll^ : solid lines), and ratio T^/Tp after one 
cooling cycle, steps 1-6, normalized by its initial value T^/Tp 
(Eq. (|22f) : dashed lines), both as a function of the fermionic 
fraction x. The colors correspond to the JILA parameters 
(red) , the MIT parameters |l3 (blue) and the best case Vhf — 
1 (black). 



does not lead to a simple compact expression. However, 
by setting T'^ = = Tiso we find the initial tempera- 
ture for which an adiabatic sweep is also isothermal (in 
the sense that the temperature endpoints are equal): 
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If > Tiso the gas is heated by the adiabatic sweep, 
while a starting temperature T* < Tiso leads to cooling. 
Assuming extremely low temperatures, such that the en- 
tropy can be taken to be linear in T, a simple relation 
can be obtained for AN > 0: 



Tf 



,2/3 



(2a;- 1)2/3 + r6/(l-a;)2/3' 



X > 0.5. (19) 



In the case of equal atomic populations, the relation be- 
comes: 



(20) 



The same expression holds in the case of AN < 0, where 
the ratio /T^ is independent of |AA^|, as long as there 
is a condensate in both initial and final states. This is 
the regime where experiments are currently operating. 

Figure El plots the dependence of the final temperature 
on the initial temperature for x = 1 /2 for the JILA trap 
parameters 15j in comparison to the low- and high-T 
limits (pn|) and (P3|) . respectively. For these parameters 



Eq. (0 give Tiso — 0.20Tf, which is seen to provide a 
good estimate of the value of T* for which = T*. 

Cooling Cycle-We plot Eq.lUHI) in Fig. El As IS ev- 
ident, jT^ depends strongly on the fermionic fraction 
X. In particular, there is an optimal fraction, Xopt = 
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2r^j], which maximizes cooling; it corre- 
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sponds to the minimum of the /T^ curve. 

Since Xopt > 1/2 the following cooling cycle is possible: 

1. Start with an atomic gas at Cres/^B^F ^ 1- Cool the 
gas to T < Tiso through standard cooling techniques. 

2. Alter the fermionic fraction to Xopt- 

3. Sweep eres adiabatically to — Cres/^B^F ^ 1, forming 
the maximum number of molecules possible. 

4. Remove the free atoms from the system. 

5. Adiabatically increase eres into the dissociated regime. 

6. Continue the cycle by repeating steps 2 to 5 until the 
desired final temperature in reached. 

Assuming T <C Tp we can use (|T^ and to calcu- 
late the final temperature of the system after n cycles: 
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where we have neglected a small decrease in temperature 
after suddenly decreasing the populations in steps 2 and 
4. Notice that values of r5/ closer to unity produce better 
cooling; this can be attained by appropriate choice of 
wavelength of the trapping laser jl^. 

As is the case for evaporation, this cooling scheme 
leads to particle loss from the system in steps 2 and 
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4. This results in a reduction of the Fermi tempera- 
ture Tp. Hence the figure of merit for increasing phase- 
space density is the decrease of the ratio of the temper- 
ature T to the Fermi temperature Tp at the end of the 
cycle. A simple calculation shows that after n cycles 
= T^[(l - The expression for T/Tp as a 

function of the number of cycles performed is then: 
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(22) 



Fig.inishows this expression for n = 1 and Fig.^shows an 
example of the change in T/Tp, as a function of the num- 
ber of cycles performed. For x > Xopt the cycle results 
in a net increase of T/Ty^ and for x = Xopt the phase- 
space density remains unchanged. We find that there is 
a fermionic fraction, Xopt < ^opt^ which maximizes the 
reduction of T/Tp; this corresponds to the minimum of 
the respective curve depicted in Fig. |S1 We plot it as a 
function of r5/ in the inset of Fig.^ 
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FIG. 4: Ratio of absolute to Fermi temperature after n cycles 
operated at The numerical results (symbols) are 

compared with the analytical estimate (|22|) . which neglects 
a small temperature decrease that occurs in steps 2 and 4 
(lines). Results are shown for the JILA parameters 15] (as- 
terisks and dashed line) and for the best case r^/ = 1 (circles 
and solid line). Inset: function of the frequency ratio 

The efficacy of the proposed cooling scheme is contin- 
gent on losses from three-body collisions being minimal. 
Presently, no estimates of their rates exist for the systems 
under consideration. Additionally, phase separation may 
become an issue, as the interspecies atomic scattering 
length diverges on resonance 18]. 

Conclusion-We have presented the phase diagram 
of an ideal mixture of bosonic and fermionic atoms in 
thermal and chemical equilibrium with heteronuclear 
fermionic molecules. A cooling cycle has been identified 



that exploits the mechanism of adiabatic atom-molecule 
interconversion. This cycle could provide a useful com- 
plement to existing cooling techniques. Its practical effi- 
ciency will be limited by effects of atom-atom and atom- 
molecule interactions, which lie outside the scope of the 
present treatment, and must be determined by experi- 
ment. 
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